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I. 


INTRODUCTION 


Problems  of  analysis  and  synthesis  of  radiating  systems  require 
appropriate  functions  in  which  to  expand  a  source  distribution  and 
associated  radiated  field.  To  be  most  useful  it  is  desirable  that  such 
functions  be  suitable  for  source  and  field  representations  simultan¬ 
eously,  that  they  be  complete  and  that  they  be  orthogonal  in  some  sense 
over  both  the  source  and  field  region  of  interest.  In  the  familiar 
case  when  the  regions  of  the  source  and  the  field  coincide  with  coordin¬ 
ate  surfaces  of  coordinate  systems  in  which  the  Helmholtz  equation  is 
separable,  the  corresponding  eigenfunction  representation  is  valid  for 
both  the  source  and  the  field.  The  characteristic  modes  of  loss-free 
bodies  of  arbitrary  shape,  introduced  by  Garbacz  [1,21,  and  elaborated 
by  Harrington  and  Kautz  [3,41  give  another  possible  representation  where 
the  source  and  field  observation  regions  coincide,  but  the  orthogonal 
properties  of  these  functions  apply  only  to  the  body  surface  and  the 
sphere  at  infinity. 

In  this  report  we  treat  the  more  general  case  where  the  region 
of  observation  dens  not.  necessarily  coincide  with  the  source  region, 
resulting  in  orthogonal  properties  over  more  general  regions  than  the 
body  surface  and  the  sphere  at  infinity.  An  investigation  of  the  validity 
of  Parseval's  relation  for  a  more  general  operator  equation  than  a  Fourier 
transform  leads  us  to  an  eigenvalue  equation  of  the  Hermitian  iterated 
operator,  whose  solutions  shall  be  called  eigenscurces  and  the  radiated 
fields  of  which  shall  be  called  eigenfields.  An  cigensource  and  an 
associated  eigenfield  each  satisfies  an  orthogonality  property  simul¬ 
taneously  but,  in  general,  over  different,  regions  of  space.  The  sets 
of  eigensources  and  eigenfields  can  then  be  used  as  convenient  basir 
functions  with  which  to  solve  various  problems  systematically. 

[xamples  are  developed  for  the  optimization  of  an  array  under 
given  constraints  and  for  the  determination  of  aperture  distributions 
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transferring  maximum  power  to  a  second  aperture.  In  the  latter  case, 
the  eigenfunctions  arising  from  the  present  theory  are  shown  to  be  the 
prolate  spheroidal  functions.  In  the  case  where  the  eigenequation  is 
developed  from  a  current  source  flowing  on  a  closed  surface  and  the 
component  of  the  corresponding  electric  field  tangential  to  the  same 
surface,  we  are  led  to  ei gencurrents  and  eigenfields  which  differ  from 
characteristic  modes  (which  involve  the  same  two  quantities);  whereas 
a  characteristic  current  and  a  characteri stic  field  display  differing 
amplitude  distributions  but  maintain  a  constant  phase  relationship  over 
the  surface,  the  ei gencurrents  and  eigenfields  introduced  here  display 
similar  amplitude  distributions  but  maintain  a  complex  cot  jugate  phase 
relationship  over  the  surface  f 5 1 .  Only  when  the  surface  corresponds 
to  a  coordinate  of  the  coordinate  systems  where  separability  applies 
do  the  character i sti c  functions  and  the  eigenfunctions  of  the  Hormitian 
iterated  operator  coalesce  to  become  one  and  the  same  such  as  in  the 
case  of  the  circular  conducting  cylinder  and  conducting  sphere.  Examples 
are  discussed  involving  the  infinite  rectangular  conducting  cylinder, 
and  the  linear  conductor. 

II.  EIGENFUNCTIONS  ASSOCIATED  WITH  SOURCE  AND  FIELD 

We  consider  a  source  s  and  a  resulting  field  f  all  in  a  linear 
medium  which  are  related  by  an  integral  operator  G  as 

f  =  Gs  .  (1) 

When  the  operator  equation  (1)  happens  to  be  a  Fourier  transform  (when 
for  example  s  is  a  planar  source  and  f  is  the  corresponding  far  field 
on  the  sphere  at  infinity)  we  know  that  f  and  s  satisfy  a  Parseval ' s 
relation,  i.e.,  the  norm  of  f  is  equal  to  the  norm  of  s.  In  the  fol¬ 
lowing  section  we  inquire  whether  or  not  this  same  result  applies  to 
the  more  general  operator  equation  (1). 
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A .  Eigenvalue  Equation  of  a  Hermitian  Iterated  Operator 


To  begin,  let  us  define  the  inner  product  of  two  sources,  (s.,s, 
and  the  inner  product  of  two  fields  (f^.f^)  as  follows: 


(Sj,S  -)  =  /  s -(kx)s^(kx)d(kx) 

JfiS  Rs  J 


(?) 


=  f 

R 


f 


w(kx)fj(kx)fj(kx)d(kx) 


(3) 


where  w(kx)  is  a  real  weight  function  (often  unity)  and  and  R^  are 
the  regions  where  the  sources  are  distributed  and  the  fields  are  ob¬ 
served,  respectively.  The  superscript  c  denotes  complex  conjugate  and 
k  =  i s  the  wavenumber  in  terms  of  the  wavelength  x. 


Using  Equation  (1),  we  obtain 


(frV, 


•  <*,.«*<*  j) 


(4) 


where  <3*  is  the  adjoint  operator  of  The  iterated  operator  is 
Hermitian  and  is  denoted  henceforth  by+ 

«  h  (5) 


+More  explicitly,  since  H  is  an  integral  operator,  Equation  (1)  is  of 
the  form 


f  ( k  x  )  =  (is  ( kx 1  ) 


f  G(kx,kx')s(kx')d(kx') 
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(Footnote  Conti nued) 


where  G  is  the  integral  kernel.  Therefore  Equation  (3)  becomes 


=i"  w(kx)jJ  G  ( k  x ,  k  x '  )s.  ( k  x '  )  d  ( k  x  '  )i 

J  fRf  L\  J 


J  GC(k:<,kx")s^(kx")d(kx") 


d  (  kx  ) 


-J  s.(kx")</  ij  w(kx)G( kx ,kx‘ )GC (kx,k x"  )•', (  kx ) 

I  In  i  fi 


R  1  R  , 
SL  f 


S  T  (  kx  11  )d(  kx  )  \d(  kx  '  ) 
J  J 


r  1 

=  /  s  .  (kx*  )\  J  i  /  wf(kx)GC(kx,kx’  )G  (kx ,  k  x"  ;d  (kx  1 !  s  ,( k  x" )  d(  kx" ) 


l 


J'Rs  'Rs~' 


b  I  J  » 

lRsLRf 


4  C 

d  ( k  x '  } 


where  we  have  used  Equation  (4)  on  the  left  side  of  the  above  equation. 
Understanding  that  the  inner  product  ovei  the  source  region  is  defined 
by  Equation  (2)  and  comparing  the  two  sides  of  this  equation,  we  see 
that 


tf3j(kx'  )=ii*;iSj(kx‘ 


-/  \I 
\lr 


w(  kx)G'"(  kx ,kx '  )G( kx  ,kx" )  d(  k>  ; 
f 


kx") 


In  this  expression  we  recognize  that  wC(kx)  =  w(kx)  because  the  weighting 
function  w(kx)  is  real.  It  is  more  usual  to  define  an  iterated  operator 
without  the  weighting  function. 


If  G  is  the  Fourier  operator,  f*  reduces  to  an  identity  operator 
and  Equation  (4)  is  Parseval's  relation.  For  a  similar  relation  to 
hold  in  our  case  with  a  more  general  operator,  the  source  must  satisfy 
an  equation  of  the  form 
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Hsi  zVi 


(6) 


whence  Equation  M)  becomes 


(f.,f 


MW, 


(/) 


We  recognize  Equation  (6)  as  an  eigenvalue  equation  of  /;  with  eigenvalues 
v ■  and  eigenfunction  solutions  s-  which  satisfy  the  orthonormal i ty  con¬ 
dition 


T  J 


=  <S 


i  J 


(9) 


where  is  the  Kronecker  delta.  From  Equations  (4),  (7)  and  (8)  we 
al so  have 


(fi’V 


R, 


(s 


V .  f-  .  . 

i  ij 


(9) 


We  shall  call  {s-j  the  ei gensources  defined  on  and  (f-|  the  eiqenfields 
defined  on  R^.  The  iterated  Hernitian  operator  property  of  p  assures 
us  that  the  l-j.j  are  real  and  positive  semidefinite  and  that  the  { s ^ } 
are  orthogonal  and  complete  on  R$  if  i-  is  not  improperly  singular. 

The  same  may  be  said  of  the  |  f  ^ } .  We  shall  order  the  eigenvalues  {v. } 
such  that  their  values  decrease  with  increasing  index  i.  This,  together 
with  Equations  (8)  and  (9),  ini|  lies  that  the  ratio  of  the  norm  of  f 
to  the  norm  of  s  is  largest  for  the  lowest  eigenvalue  and  equal 
to  it  numerically. 


Now  that  we  have  defined  pertinent,  quantities  in  qenpral  terms, 
let  us  consider  a  few  supple  hut  useful  examples  of  the  above  theory. 

+  V'ilt 

In  the  following  sections  an  e  v  time  convention  is  assumed  and 
suppressed. 


B .  Eigensources  and  E  i  genpatterns  of  Arrays 

Figure  1  shows  an  N-olement  linear  array  of  isotropic  point  sources 
with  uniform  spacing  d.  We  consider  the  source  s  to  be  the  N  array 
elements  with  excitation  strengths  s-j,  s0  ...  defined  on  the  array 
proper,  and  the  field  f  to  be  the  corresponding  far-r3diated  field  pattern 
function  f(9)  defined  on  the  sphere  at  infinity.  In  this  case  the 
operator  Equation  'll  becomes 


f(G)  =  7  GJn)s  . 


i  0  n  ’ 
n=  1 


(n\  .  jnkdcosr 


n  -  1 ,  ? . . .  N 


and  Equations  (?)  and  (3)  become 


,s,s  „  =  s  s 

’  'R  s  n  n 
s  n  - 1 


(f,f)R  =  2m  /  f(0)fC(e)si.i0d0  (13) 

f  b 

The  weight  function  w  in  this  case  is  2u$in-T.  The  eigenvalue  equation 
(6)  reduces  to  the  algebraic  equation 


H(n,m)s-  -  \>  s.  ,  i,n  =  l,?,...N 

v  ’  '  lm  i  in 


where  the  operator  a  reduces  to  an  NxN  matrix  H  with  elements 


H( n.ro) 


2v  f  G0(n)G^(m)s  iivhin 
o 


.  s  i  n  i  n-m  4.d 

T  71  - 7  V  -I.:- 

(  n  -  in  ( k  d 


n,m= 1 . . N 
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,iii4i,j,yL 


The  eigenpattern  f .(0)  corresponding  to  the  ith  eigensource  s.  is 

fi'°>  =  5,  Go(n,sin  ■  06) 

The  sets  { s  • }  and  { f  }  are  useful  in  at.  least  two  ways.  First, 
the  eigensource  with  the  largest  eigenvalue  is  the  source  which  radiates 
the  largest  power  while  the  norm  of  the  source  is  fixed,  i.e.,  the  array 
has  the  highest  radiation  efficiency.  Secondly,  we  can  synthesize  a 
given  arbitrary  pattern  f(0)  in  a  least  mean  square  sense  by  choosing 
the  N  source  excitations  in  Equation  (11}  t'  he 


N 

y 

i  ■  1 


(M-T) 


i  n 


nl,\. .  .N. 


(17) 


Fi ymes  ?  and  3  sliow  calculated  eifjensojrr.es  and  eigenpattern 
magnitudes  for  N-S,  r!D.4x. 


C .  Focussing  the  Field  of  an  Array 

As  a  second  example  consider  the  same  linear  array  discussed  above, 
but  whose  field  is  observed  at.  a  point  P  in  the  near  field  region  as 
pictured  in  Figure  4.  Here,  the  appropriate'  operator  reduces  to  a 
1 xN  matrix  with  elements 


-  J'<r 


G ,  ■  n  )  -  r  — 

I  '  M' 


n - 1 N 


(  IB) 


where  the  r  are  def  ine.l  in  F  igure  4,  (..id  the  operate'  /i  reduces  to 
an  NxN  matrix  witli  eleineiit1. 


H ( n.i'i) 


j  M  r 


n 


r 


n 


r* 

in 


(i  ,m  1  ,  ?  .  .  .  N 


(  19) 


n 


The  norm  of  s  remains  the  same  as  Equation  (12)  but  because  region 
in  this  case  is  the  near  field  point  P  the  norm  on  f  becomes  simply 

(f.f)p  ^  f(P)fC(P)  •  (?0) 

It  i s  obvious  that  the  eiqensource  with  the  largest  eigenvalue 
yields  the  largest  field  at  point  P  while  the  norm  of  the  source  is 
fixed.  Figure  5  shows  this  case  of  an  optimally  focussed  1  inear  array 
with  N  =  5,  d~0.5X  and  P  located  at  (2X,1.5X)  --  (xp,Zp).  Wo  comment  that 
the  optimum  field  at  ?  is  larger  by  0.11  dB  than  is  the  field  obtained 
there  from  an  array  et  co-phascd  elements  with  uniform  amplitudes. 

D.  Eiqensources  and  Eigenfields  of  a  Slit  Aperture 

As  an  example  of  a  continuously  distributed  radiating  system, 
consider  two  identical  planar  parallel  slits  of  width  ?a  and  separation 
b,  as  Shown  in  Figure  6.  The  left  slit  is  considered  to  form  a  source 
and  the  right  slit  serves  to  observe  the  field  due  to  the  source.  For 
such  a  two  dimensional  problem  (i  is  the  integral  operator  !  h'| 

ka 

f ( k  x  p)  =  /  G(kx^,kx^)s(kxj)d(kx|)  (?1) 

L  -ka 

with  kernel 

i  J!  ( x  ] _x?) 

J  A  rr  -  jk — — •  - 

G(kx?,kx,)  -  'J-  e  -  (?.?) 


The  Fresnel  approx  imat  i on  has  been  assumed  in  Equation  (?2).  Thp 
Hermitian  iterated  operator  a  is  an  integral  operator  with  kernel 


Z_Sn 


OPTIMUM 

UNIFORM 
AND  COPHASE 


Figure  5.  Optimum  source  distribution  nf  five  element  linear 
array  antenna  focussed  in  tlic  near  field  region  compared 
with  uniform  amplitude  and  cophase  distribution. 
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ka 


H(kx^kx'j)  -  /  GC(kx2,kx^)  G(kx2,kx".  )  d(kx.,) 
—  k  d 


d 

7Tb 


hjk 


x'2_x'J?  1 
X1  X1 

— 2b " 


sin  ^(x)-x'p 


(23) 


Defining  a  modified  source  function  r(kx-|)  and  a  real  symmetric  kernel 
h(kx^,kx'j)  by 

(kx^2 

■j  2ktT 


r(kx^ )  =  e 


s(kx1 ) 


(24) 


h(kx).kx'i)  h  fb 


sin  |(kxj-kxp 
|(kx|-kxp 


(25) 


the  eigenvalue  equation  (6)  becomes 


ka 


j”  h(kx  ^  ,k  x'j  )r  •  ( kx'|  )d(  kx'j )  =  ^r.fkxj)  ,  i  =  0,1,2,...  (26) 

-ka 

Once  this  equation  is  solved  for  {r  -  } ,  the  set  of  eigensources  {s  ^ } 
and  eigenfields  { f  • }  are  given,  respectively,  by 


s  .  ( k  x ,  1  =  N  •  r  .  ( k  x  n )  e 
i  i  i  i  r 


(kx,) 


(27) 


and 


ka 


f^kx,,)  =  Ni  J  G(kx?,kx1)  s^kx^)  d(kx^) 
-  k  a 


-J, 


=  N  ■ 


( k  Xp 

~?Fb 


\2 


-A  ! 


ka 


k  x  ^  k  x  i 
+  kb 


-  ka 


n ^ ( k x i )  d(kx  1 ) 


(28) 

where  N.  are  convenient  normalizing  constants.  With  a  change  of 
vari ablos 


ka 

T 

(29) 

a 

b 

=  a 

(30) 

kx} 

=  t 

(31) 

kx" 

:  S  , 

(32) 

equation  (26)  may  be  rewritten  in  the  form 


T  /? 

f  sinQ(t-s) 

_'T /2  77  (T^sJ 


r  .  ( s )  ds 


’iri (t) ,  i  =  0,1,2... 


(33) 


This  equation  has  been  studied  extensively  by  Slepian  and  Poliak  <  7'! 
who  show  its  solutions  to  be  the  angular  prolate  spheroidal  functions, 
r  i  ( t )  =  S  j(c,t)  and  its  eigenvalues  to  be 


v.j(c)  - 


2c 


,d), 


01 


(c.,1) 


1  =  0,1,2, 


(34) 
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where  rH)(c,1)  are  radial  prolate  spheroidal  functions.  The  parameter 

01  9 

O 

c  =  ~  or>  according  to  Equations  (29)  and  (30),  c  =  —  ■  =  2rF  where 

F  is  the  Fresnel  number 

F  =  —  .  (35) 


Asymptotic  expressions  for  A^(c)  for  small  and  1 arge  values  of  c  (i.e., 
small  and  large  values  of  the  Fresnel  number)  with  index  i  being  fixed 
have  been  derived  by  Slepian  and  Sonnenblack  [3). 

For  the  normalization  condition  expressed  by  Equation  (8)  to  hold, 
the  ith  eigensource  is 


s  .  ( k x , )  =  v',,/2(c)  e 

1  1  i 


{ Is  x  -J  )  £ 


Soi<c’kxl) 


and  the  corresponding  eigenfield  given  by  Equation  (28)  is 


fi(kx2)  =  v-^(c)  e 


,PV2  ™  nn 


f  ( k  Xp ) c 

.[/•ihB--i-rJSo((Cilt  Js  4m(c,„].,37) 

Since  the  fS  -}  form  a  complete  orthogonal  set  on  either  slit,  a  source 
and  its  field  are  conveniently  expanded  in  terms  of  (s.)  and  { f  ^ } ,  re¬ 
spectively.  The  eigenvalue  v.(c)  corresponds  to  the  ratio  of  the  power 
received  to  the  power  transmitted  in  the  ith  mode.  Furthermore, 
the  i=0  mode  has  the  largest  eigenvalue,  vQ(c),  from  which  we  conclude 
that  power  transfer  between  apertures  is  maximized  in  this  lowest  order 
mode.  Using  this  fact  and  the  fact  (8)  that  for  small  c  (i.e.,  small 

n. 


ti  nn 


v0(O  =  4F, 


i 


we  can  state  that  for  two  distantly  separated  slits  the  received  power 
cannot  exceed  4F  times  the  transmitted  power.  For  large  c  (i.e.,  large 
F)  it  is  known  T 8)  that 

vQ(c)  *  1  -  4  ^fc  e'2c  -  1  -  4n  /ZF  e'4T,F  (39) 

which  has  a  value  only  slightly  less  than  1  for  F  values  larger  than 
F  =  1/2. 

It  is  worthwhile  to  note  that  the  lowest  order  (optimum)  source 
has  a  positive-going  parabolic  phase  disti  sbution  (Equation  (27))  which 
tends  to  compensate  for  the  diverging  effect  of  a  finite  slit;  the  field 
achieved  by  this  optimum  source  has  a  negative-going  phase 
distribution  (Equation  (28)). 

III.  EIGENVALUE  EQUATION  WITH  HERMITIAN 

WEIGHT  OPERATOR 

It  has  been  shown  in  the  previous  sections  that  the  eigenvalue 
equation  of  the  Herinitian  iterated  operator  yields  eigensources  anti 
eigenfields  which  can  form  useful  basic  sets  in  which  to  expand  arbitrary 
source  and  field  distributions.  It  was  also  shown  that  an  optimum  source 
can  be  determined  which  produces  maximum-field  over  the  designated  field 
region  under  the  condition  that  the  norm  of  the  source  is  fixed.  In 
this  section  we  extend  the  theory  to  encompass  more  general  conditions 
by  replacing  Equation  (6)  by  an  eigenvalue  equation  of  the  Hermitian 
iterated  operator  which  includes  a  Hermitian  weight  operator,  i.e., 

i is  ■  =  V.  ?is  ■  .  (40) 

Here  we  assume  that  the  weight  operator  is  Hermitian  and  positive 
definite  in  which  case  it  can  be  shown  that,  the  solutions  of  Equation 
(40)  satisfy  the  relations 


1 


■  5io 
's 


V,Sj>  =  Vij 

s 


(41) 


(4?) 


Equations  (40),  (41),  and  (42)  are  extensions  of  Equations  (6),  (8), 
and  (9),  respectively.  Thus,  we  can  interpret  the  eigensource  with 
the  largest  eigenvalue  as  being  the  optimum  source  in  the  sense  that 

(s,h$)r  is  largest,  (s,hs)r  being  fixed  (rather  than  (s,s)R  being 
s  s  s 

fixed  as  developed  in  Section  II). 

In  work  to  follow  on  arrays  it  will  be  useful  to  articulate  the 
operators  in  Equation  (40)  in  the  form  of  matrices.  If  U  is  the  orthog¬ 
onal  matrix  which  diagonalizes  H  and  p  is  the  diagonalized  matrix,  i.e., 

U*H  U  =  p  = 

then  Equation  (40)  transforms  to 

(U'1/2  U*HUP‘1/2)(p1/2U*si)  =  v.(p1/2U*Si),  (44) 

or 


(43) 


=  v; 


whose  eigenvalues  satisfy  the  orthogonal i ty  condition 


(45) 


IQ 


■5. 

ij 


(46) 
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It  is  easy  to  verify  that  Equations  (45)  and  (46)  are  consistent  with 
Equations  (42)  and  (41),  respectively. 

A 

A  second  useful  case  arises  when  H,  in  addition  to  being  Hermitian 
and  positive  definite,  is  also  factorable  as 

H  =  A  *A  .  (47) 


Equation  (40)  then  transforms  to 


-1* 


)  (A  s  i )  « 


V' 


(40) 


or 


l  1  1 


(49) 


whose  eigensolut ions  satisfy  the  orthonormal i ty  condition 


(S'l,  S'!)  =  A 


1’  J 


R 


1 J 


(50) 


It  is  easy  to  verify  that  Equations  (49)  and  (60)  lead  back  to  Equations 
(42)  ana  (41),  respectively.  In  addition  it  ran  be  shown  '  9’,  that  when 
li  is  factorable  the  minimum  eigenvalue  is  (in  matrix  notation) 


v  ■  =  (  A  H  AM 

min' 


(51) 


and  the  corresponding  eiqensource  is 


S.  - 
min 


(52) 
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A .  Maximizing  Eridfire  Radiation  of  an  Array 
While  Fixing  the  Radiated  Power  in  a 
Specified  Angular  Sector 

A  linear  array  of  N  isotropic  point  sources  equispaced  a  distance 
d  along  the  z-axis  as  sketched  in  Figure  1  radiates  into  the  conical 
sector  about  endfire  shown  in  Figure  7.  The  power  radiated  is  given 

by 


P  =  2" 


m=  1 


N 

y 

n=l 


Vn 


i 

e. 


i'  (m-n)kdcosO 


sinO  dO 


where 


N 

c 

m=l 


N 

V 

n- 1 


H(m,n) 


(53) 


H(m,n)  =  2n  e 


-j(m-n) 


kd 


(cosC‘0-l) 


s  inkm-n) 


kd 

7~ 


( l  +  cos0o)j 


i  \  kd 
(m-n)  — 


(54) 


H(m,n)  are  the  elements  of  an  NxN  Hennitian  matrix  A.  It  is  clear  that 

A 

P  appears  in  the  form 


P  =  (s,Hs)r  (55) 

s 

which,  if  we  constrain  P  to  be  un;ty,  corresponds  to  Equation  (41). 
Similarly,  the  power  density  radiated  in  the  endfire  direction,  G=r,  is 


5  =  ( s , Hs )  R 

s 


(56) 


?! 


where  H  is  an  NxN  matrix  whose  elements  are 


H(m,n)  =  e"J'(m'n)kd  .  (57) 

The  optimum  source,  that  is,  the  source  winch  maximizes  the  end- 
fire  power  density  S  while  constraining  the  radiated  power  P  to  a  specified 
sector  is  the  eigensource  with  the  largest  eigenvalue  satisfying  Equation 
(40)  with  H  and  H  specified  by  matrix  elements  given  in  Equations  (54) 
and  (57),  respectively.  Field  patterns  for  an  example  with  N=5,  d=0.3X, 
Oor30°(30°) 120°  are  shown  in  Figures  8-10. 

IV.  EIGENVALUE  EQUATION  OF  PERFECTLY  CONDUCTING  BODIES 

Garbacz  '  rvort  .gat.  d  ,■  oener.l i/ed  exp-vn  ion  tor  the  field 
radiated  or  scattered  by  a  loss-free  obstacle  of  arbitrary  shape  S  in 
terms  of  character i st i c  functions  defined  in  Ref.  (!’.  He  nnt^d  for  the 
perfectly  conducting  obstacle  that  each  characteristic  function  is  as¬ 
sociated  with  a  real  character istic  current  defined  on  S  which  gives 
rise  to  equiphase  fields  throughout  the  volume  occupied  by  the  obstacle. 
This  observation,  though  not  proven  in  a  general  way,  led  him  to  obtain 
a  character i st i c  equation 


-  Ai 

(58a) 

-  J,  =  (1+jV  k  ^ 

(58b) 

where  J-  is  the  ith  characteristic  current  associated  with  the  ith 
character i Stic  value  \]  and  R  and  X  are  the  Real  and  Imaginary  Hermitian 
parts 


(59) 


X 


(60) 


a 


o 


F igure 
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8.  Optimum  field  radiation  patterns  of  an  array  whose  radiated 
power  in  a  sector  of  angle  (o  ,n)  is  minimized  while  the 
endfire  radiation  intensity  is  fixed. 


of  the  impedance  operator  z  =  R+jX  of  the  obstacle  !  ldl.  The  function 
Ei  is  the  ith  characteristic  field  component  tangential  to  the  surface 
S.  The  characteristic  currents  (J.)  simultaneously  yield  orthonormal 
characteristic  radiation  patterns  while  they  themselves  have  the  ortho¬ 
normality  properties 

(RJ.,J.)  =6  (61) 

J  R 

s 

-xiau  <62) 


where  the  inner  product  is  still  defined  by  Equation  (?). 

We  next  proceed  to  discuss  the  set  of  eigenfunctions  of  the 
Hermitian  iterated  operator  as  a  set  suitable  for  expanding  fields 
radiated  or  scattered  by  a  perfectly  conducting  obstacle.  To  begin, 
we  note  that  the  impedance  operator  :  relating  current  on  5  to  tangential 
field  on  S  is  an  example  of  the  operator  6  described  in  Section  II  and 
so  we  can  define  eigensources  and  eigenfields  in  the  manner  already 
described.  In  this  case,  the  eigensources  are  currents  K .  on  the  obstacle 
surface  and  the  associated  eigenfields  are  components  of  the  electric 
fields  tangential  to  the  same  surface  (making  Z  unbounded).  The  eigen- 
currents  are  solutions  of  Equation  (6) 


k,  -  v.  Ki 


(63) 


satisfying  the  orthononnal ity  conditions,  (Equations  (8)  and  (9)) 


“VKj>  ‘sij 

s 


(64) 


(Ki.-*’Kj)r  =  V 


1 J 


(65) 
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Note  that  the  character! Stic  values  X ■  are  always  real  and  can  take 
values  -  oo  £  while  the  eigenvalues  v  ■  are  always  real  and  can  take 

values  0  <  u-.  The  characteristic  values  X^  =  -  ®  correspond  to  charac¬ 
teristic  currents  which  resonate  the  interior  volume  of  the  obstacle 
while  eigenvalues  =  0  correspond  to  this  condition. 

When  the  ambient  medium  is  reciprocal,  i.e.,  7*-Z  ,  Equation  (63) 
can  be  reduced  to  a  simpler  and  instructive  form.  Operating  on  Equation 
(63)  with  7, 


z:*(z  k.)  =  Vi(:  Ki)  , 

:zc(zK^  =  vi(z|<i)  (66) 


and  comparing  this  equation  with  Equation  (63),  conjugated. 


ZZC  =  v.  (67) 

we  see  that  z  K;  and  K1^  satisfy  identical  eigenvalue  equations.  This 
being  so,  they  must  be  related  by  a  multiplicative  constant  c^,  i.e., 


^i  =  Ci  S 


E!  =  Ci  *1 


(68) 


where  E-  is  the  ith  eigenfield  component  tangential  to  the  obstacle 
surface.  The  magnitudes  of  c-  are  fixed  by  Equations  (64)  and  (65) 
to  be 


(69) 


while  the  phases  of  Cj  are  not  unique,  being  related  to  the  arbitrary 
phases  of  K.  by  Equation  (60).  Thus,  we  can  choose  the  phases  of  c. 
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1/2 

for  convenience.  One  possible  choice  is  c.  =  jv/  ,  purely  imaginary 
quantities,  in  which  case  Equation  (68)  can  be  rewritten 


where  K1.  and  K‘.‘  are  the  real  and  imaginary  parts  of  K^,  respectively. 

1/2  ^Yi 

Another  possible  choice  is  c.  =  v-  e  where  are  chosen  so  that 
the  phase  of  is  the  negative  cf  the  phase  of  K . ,  that  is,  the  eigen- 
current  and  assorted  eigenfield  £  ■  are  scaled  conjugate  functions 
over  S.  This  latter  choice  for  c^  is  general  and  provides  an  interesting 
contrast  between  eigenmodes  of  the  Hermitian  iterated  operator  and 
characteristic  modes;  namely,  the  real  characteri sti c  current  J.  and 
its  associated  tangential  field  £.  maintain  a  constant  phase  difference 
over  the  surface  S  while  their  amplitude  distributions  differ  in  general 
over  S;  on  the  other  hand  the  generally  complex  eigencurrent  and 
its  associated  tangential  field  maintain  the  same  (scaled)  amplitude 
distributions  over  $,  while  their  phases  vary  over  S  in  such  a  way  as 
to  make  them  conjugate  functions. 

For  those  special  conducting  shapes  S,  such  as  the  sphere,  the 
circular  cylinder,  the  elliptic  cylinder,  etc,  which  correspond  to 
constant  coordinate  surfaces  of  the  separable  coordinate  systems,  the 
characteristic  currents  J-  and  ei gencurrents  K.  become  identical  (except 
for  different  real  normal izat ion  constants)  and  their  associated  fields 
E-  and  become  identical  (except  for  different  real  normalization 
constants).  Thus,  the  modal  currents  and  fields  on  S  track  in  both 
amplitude  and  phase.  In  such  special  cases,  the  ratio  of  the  ith  modal 
field  and  its  current  becomes  a  complex  (in  general)  number  which  remains 
fixed  over  the  entire  surface  and  may  be  interpreted  as  the  ith  modal 
impedance  Z-  =  8^  +  jX..  Then  Equations  (58),  (68)  and  (69),  together 

1/?  JYj 

with  the  choice  c}  -  v.  “  e  lead  to 
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kW,  e 


X  . 

j  tan  (p^-) 


or 


(R2+x?) 

i  i 


and 


Ai 


Xi 

r 


(71) 


(72a) 


(72b) 


which  indicate  the  connection  between  the  eigenvalues  v  of  the 
Hermitian  iterated  operator  and  the  character i st i c  values  X-, 


.  ,  -1. 

-v  J  tan  X, 

ci  =  RiVl  +  X:  e 


or 


?  ? 

v.  =  R.  (1  +  x "  1 

i  i  v  i 


(73a) 


(73b) 


As  examples,  Equation  (70)  has  been  applied  to  a  thin  straight 
wire  and  infinite  cylinder  of  rectangular  cross  section,  using  the  method 
of  moments  ;  10 1  to  arrive  at  matrix  representations  for  p  and  X.  Values 
of  v-  are  shown  in  figures  11-13  as  functions  of  electrical  size  of 
tiie  obstacle.  As  expected  in  the  case  of  the  cylinder,  where  we  pass 
through  resonances  of  the  interior  region  we  notice  that  certain  v  • 
become  zero  for  kb  corresponding  to  cut-off  frequencies  of  associated 
waveguide  modes.  Those  frequencies  agree  with  theoretical  values  within 

0.5%  for  the  TM  modes  and  within  0.24%  for  the  TE  ""'les,  one  quarter 
of  the  cylindrical  periphery  having  been  divided  into  IT  parts  for 
application  of  the  method  of  moments. 
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Figure  12.  Eigenvalues  of  TM  modes  for  a  rectangular  cylinder  with 
side-lengths  b  and  4b/9.  The  numbe1"  of  subdivisions  in  the 
computation  by  the  moment  method  is  13  for  one  quarter 
of  the  periphery. 
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Figure  13.  Eigenvalues  of  TE  modes  for  a  rectangular  cylinder  with 
side-lengths  b  and  4b/9.  The  number  of  subdivisions  in  the 
computation  by  the  moment  method  is  13  for  one  quarter 
of  the  periphery. 
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V.  CONCLUSIONS 


i 


By  extending  Parseval's  relation  to  any  combination  of  a  source 
and  a  field,  we  have  arrived  at  an  eigenvalue  equation  of  a  Hermitian 
operator  formed  by  iterating  any  general  operator  relating  a  source 
and  a  field.  The  eigensources  and  corresponding  eigenfields  whicn  are 
solutions  to  this  eigenvalue  equation  are  complete  and  orthogonal  over 
the  source  and  field  regions,  respectively.  Any  eigenvalue  of  the 
equation  is  the  ratio  of  the  norm  of  the  corresponding  eigenfield  and 
the  norm  of  the  corresponding  eigensource.  Except  when  the  operator 
is  improperly  singular,  the  completeness  and  orthogonality  properties 
of  the  o i gensol ut i ons  make  them  attractive  for  expanding  arbitrary  fields 
and  sources  associated  with  discrete  or  continuous  radiating  or  scat¬ 
tering  systems. 

A  few  examples  show  the  application  of  the  theory  to  arrays  and  to 
two  planar  apertures.  Out  of  the  latter  work  comes  the  interesting 
observation  that  the  theory  applied  to  the  operator  equation  known  as  the 
Fourier  transform  leads  to  eigenfunctions  which  are  the  prolate  spheroidal 
functions. 

When  the  theory  is  applied  to  tne  complex  impedance  matrix  of 
a  conducting  scattering  obstacle,  it  is  shown  that  the  eigencurrents 
and  eigenfields  introduced  here  are  conjugates  of  each  other,  which 
contrast  with  characteri stic  currents  and  fields  of  Garbacz,  Harrington 
and  Mautz.  Only  in  the  case  of  obstacles  corresponding  to  coordinates 
of  separable  coordinate  systems  do  the  two  modal  types  coincide. 

The  eigenfunction  tneory  is  applied  to  the  thin  linear  wire  and 
the  rectangular  cylinder  and  the  correspond i ng  sets  of  eigenmodes  are 
presented  here.  Results  encourage  the  application  of  the  theory  to 
other  shapes  and  invite  an  investigation  of  the  ei gencurrent  distributions 
as  well  as  the  eigenvalues. 
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